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Abstract 

We present continuum and lattice calculations for elastic scattering between a fermion and a 
bound dimer in the shallow binding limit. For the continuum calculation we use the Skorniakov- 
Ter-Martirosian (STM) integral equation to determine the scattering length and effective range 
parameter to high precision. For the lattice calculation we use the finite-volume method of 
Liischer. We take into account topological finite-volume corrections to the dimer binding energy 
which depend on the momentum of the dimer. After subtracting these effects, we find from 
the lattice calculation naf<i = 1.174(9) and K-rfd = —0.029(13). These results agree well with the 
continuum values nafd = 1.17907(1) and nrfd = —0.0383(3) obtained from the STM equation. We 
discuss applications to cold atomic Fermi gases, deuteron-neutron scattering in the spin-quartet 
channel, and lattice calculations of scattering for nuclei and hadronic molecules at finite volume. 

PACS numbers: 21.60.De, 25.40.Dn, 12.38. Gc, 03.65.Ge 
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I. INTRODUCTION 



In this paper, we perform the first benchmark of finite-volume lattice methods for the 
low-energy scattering of composite objects. Our results will have immediate applications 
to lattice studies of elastic neutron-nucleus scattering. In the analysis presented here we 
consider scattering between a fermion and a bound dimer composed of two fermions. In 
order to test the precision of our lattice calculations, we also repeat the same calculations 
using the Skorniakov-Ter-Martirosian (STM) integral equation Along the way we also 
provide the most accurate calculation to date for the fermion-dimer effective range param- 
eter. Some of the results presented were summarized in a letter publication [2j, and we 
present the full details of the calculations here. 

We consider two component fermions. We will refer to the two fermion components as 
spin up and spin down and consider the case when the masses are equal, = m^. We 
assume finite-range attractive interactions and consider the universal shallow binding limit. 
If R is the range of the interactions and k is the binding momentum of the dimer, then the 
shallow binding limit corresponds to kR — > 0. 

We note that much of the literature on universal three-body systems has focused on 
the Efimov effect [3| for three bosons, three-component fermions, unequal mass fermions, 
or mixed Bose- Fermi systems jJS]- For equal mass two-component fermions, however, 
there are no short-distance three-body instabilities such as the Thomas collapse [9]. Hence 
there are no relevant momentum scales other than the dimer binding momentum n, and 
all low-energy scattering parameters can be expressed as dimensionless constants times the 
corresponding power of k. In the shallow binding limit k is the same as the reciprocal of 
the fermion-fermion scattering length. 

There have been numerous calculations of the fermion-dimer scattering length. The first 
goes back to the early work of Skorniakov and Ter-Martirosian who found najd rj 1.2 [lj]. 
An overview of the Skorniakov-Ter-Martirosian (STM) integral equation method will be 
presented in our discussion below. The same value na/d ~ 1.2 has been confirmed several 



times using integral equations [lOj, lUj. A value of Kafd ~ 1.11 was obtained using an epsilon 
expansion in dimensions 12|. A recent correlated Gaussian expansion calculation obtained 
KOfd ~ 1.18(1) [13j. This agrees with integral equation studies which found Kdfd ~ 1-18 
llTll and na fd « 1.1790662349 El. 



The fermion-dimer results at shallow binding should approximately describe neutron- 
deuteron scattering in the spin-quartet channel. Experimental measurements find a quartet 
scattering length A a n d = 6.35(2) fm [l7[. This corresponds with n A a nd = 1.47(1). The 
agreement is better when expressed as fraction of the spin-triplet proton-neutron scattering 
length, 4 a n rf/ 3 a np = 1.17(1). The 30% difference between the two values gives an indication 
of higher order effective range effects. A more detailed calculation including interaction 



range effects obtains a n d = 6.33(10) fm |18|, [19J, in full agreement with experimental values. 



In contrast with the scattering length, there is only one previously reported determination 
of the fermion-dimer effective range parameter. The correlated Gaussian expansion calcula- 
tion in Ref. [13| found Krf d ~ 0.08(1). Neutron-deuteron scattering data also favors a small 
value for KTfd- However the sign of nrjd has remained an open question. In our continuum 
calculations presented here, we use the STM equation to calculate the scattering length and 
effective range parameter. We find the values Kafd = 1.17907(1) and KTfd = —0.0383(3). 

Our main focus though is to benchmark lattice calculations of the fermion-dimer scat- 
tering length and effective range parameter. For the lattice calculation we apply the finite- 



2 



volume phase-shift analysis of Liischer. We show that finite-volume topological corrections 
to the dimer binding energy must be considered in order to obtain accurate results. Once 
these topological corrections are included in the finite-volume analysis, we find that the 
lattice and continuum calculations are in full agreement. We use two different lattice 
Hamiltonian formulations, that agree in the continuum limit. As a final result we find: 

Ka }d = 1.174(9), nr fd = -0.029(13). (1) 

II. NOTATION AND FORMALISM 

Few-body systems of two-component fermions with short-range interactions and large 
scattering lengths in comparison to interparticle distances show universal properties. 
Physics in such systems does not depend on the structure of the interactions at short dis- 
tances. The problem of three two-component fermions at low energies can be described by 
a local quantum field theory whose only interaction term is a two-body contact interaction. 
In the following, we will always consider equal mass fermions with mass = = m. 
The extension to unequal masses is straightforward. The free non-relativist ic effective 
Hamiltonian in momentum space can be written as 




where a« and a] are annihilation and creation operators. In position space these operators 
can be expressed as 

<?) = (2^8 / d 3 pe^%(p), a\{r) = * J e~^a\{p). (3) 

Combining Eqs. fl2]) and ([3]) the free Hamiltonian in configuration space is then given by 




Now we introduce an interaction between the fermions via the potential 

V(r, r ') = \^2 f d 3 r f d 3 r' : a\{r)a t {f)V{r - f>J(f >j(r') : , (5) 
i,j=U J J 

where : . . . : denotes normal ordering. At low energies, the potential (JSJ) can be replaced by 
a delta-function interaction 

V(f-f') = C 5 (3) (f-f'), (6) 
and the the lowest order effective Hamiltonian for two-component fermions is 
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Co denotes the two-body coupling constant and is directly related to the fermion-fermion 
scattering length. It is assumed to be negative so that the interaction is attractive. The 
exact value of Co depends on the scheme used to regulate the short distance behavior. 

In the next step we consider the two-body and three-body systems of two-component 
fermions in Hamiltonian lattice formalism and use the Lanczos method 20] to find the 
lowest eig envalues. Further details of the Hamiltonian lattice formulation can be found 
in [il -23]. Let n denote spatial lattice points on a three-dimensional L x L x L periodic 
cube. We use lattice units where physical quantities are multiplied by powers of the spatial 
lattice spacing a\ a tt to make the combination dimensionless. The two-component fermions 
are labelled as spin-up and spin-down, the lattice annihilation operators are written as 
Of(n) and a±(n). The free non-relativistic Hamiltonian of two-component fermions with 
only short-range interaction corresponding to the Hamiltonian (T4|) on the three dimensional 
lattice is 

Ho = m^-< a i^ Qi ^ _ 9^ E E [4(™)°»(** + #) + al(n)a;(n -/})], (8) 



2m 



=1,2,3 *M=t4 



where fi is the spatial lattice unit vector. We define the spin-density operators 

Pf (n) = aL(n)a^(n) (9) 
Pi(n) = a|(n)a ; (n), (10) 

and consider two different kinds of Hamiltonians. In the first Hamiltonian we have only a 
single-site interaction. This Hamiltonian is 



H x = F + C 1 ^p t (n)p^(n) 



(11) 



We consider a second Hamiltonian using a contact interaction as well as nearest-neighbour 
interaction terms in order to eliminate the two-body effective range parameter, 



H 2 = H + C 2 ^2 /°r(")P4.(*0 

n 



£=1,2,3 n 



(12) 



latt ' 



The finite lattice spacing error in these two Hamiltonians is of order a 

The next step is to determine the interaction coefficients, C*i, C 2 and C 2 using Liischer's 
formula [24h26|. At present Liischer's formula is a standard tool in lattice quantum chro- 
modynamics and in lattice effective field theory. It relates the two-body energy levels in a 
finite volume to the S-wave phase shift. 



p cot 5 Q (p) 



7lL 



V 



Lpy 

2ir) ' 



(13) 



where S(rj) is the three-dimensional zeta-function, L is the length of the box and p is the 
center-of-mass momentum. The zeta-function in three dimensions is defined as 



6 (A 2 - ff 



4ttA 



rr 



7} 



(14) 
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For \r)\ < 1 we can expand S(rj) in powers of rj, 

S{v) = — + hm V - 4ttA 

77 a^oo L z — ' n £ — ri 



where 



ri a^oo L * — ' n z — r? 
-- + S + S x r) + S 2V 2 + S 3V 3 + ■ ■ • , (15) 



S = hm \Y d{A \ ^ - 4vrAl , S l = y — L_ (16) 



The hrst few coefficients are 



S = -8.913631, Sx = 16.532288, S 2 = 8.401924, S 3 = 6.945808, 

^ 4 = 6.426119, S 5 = 6.202149, S 6 = 6.098184, S 7 = 6.048263. (17) 

Liischer's formula does not include the contribution from higher partial waves but at asymp- 
totically small momenta we can neglect such corrections. For small momenta we have the 
effective range expansion, 

pcot<5oG»)~ — - + /o P 2 H , (18) 

(X s Z, 

where a s is the scattering length and r is the effective range. Thus the effective range 
parameters can be extracted from the finite volume energy levels using Eq. (flBl) . In terms 
of 77, the energy of the dimer is 

E = P - = ^L( 2 4-Y. (19) 



m m\L 

The interaction coefficient Cx is tuned to construct two-body binding states (dimers) 
comprised of one spin-up and one spin-down fermion of energies —1.5 MeV, —2.0 MeV, 
-2.5 MeV, -3.0 MeV, -3.5 MeV and -4.0 MeV in the large volume (L = 80) using the 
Lanczos method. In such a large volume the finite volume corrections to the dimer binding 
energy are negligible. In our calculation we take m = 939 MeV and aj^ tt = 100 MeV. To 
find the interaction coefficients of lattice Hamiltonian H 2 we proceed as follows. Setting the 
effective range to zero requires that the following relation should be satisfied near threshold 

-^rS(ri) = pcot 6 (p) ^ -- + 0(p 4 ). (20) 

The interaction coefficients C 2 and C' 2 are tuned in order to give the binding energies listed 
above for the ground state in a large volume (L = 80) and to fullfil Eq. (1201) for the first 
excited state. The plot for pcot 5o(p) versus p 2 for the first excited state is shown in Fig. CD 
The different values of p are generated by calculating for different box sizes. We note that 
pcotSo(p) has zero slope near threshold since we set the effective range to zero. Both 
Hamiltonians reproduce the same continuum limit of fermions with attractive zero-range 
interactions. The corresponding values for the interaction coefficients are summarized in 
Table E 
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FIG. 1: Plot of pcot5o(p) versus p 2 for the lattice Hamiltonian H^- 



♦, ▲, +, x , represent 



data points corresponding to the bound states of energies —1.5 MeV, —2.0 MeV, —2.5 MeV, 
—3.0 MeV, —3.5 MeV and —4.0 MeV, respectively. * denotes pcot<5o(p) at the zero mo- 
mentum limit. 



TABLE I: The values of interaction coefficients for the six considered dimers. All quantities are 
given in units of the lattice spacing a\ at t = (100 MeV) -1 . 



E d 


mC\ 


mC 2 


mC 2 


-0.015 


-4.51091 


-4.34554 


-0.30082 


-0.020 


-4.61299 


-4.45733 


-0.34273 


-0.025 


-4.70675 


-4.55883 


-0.34658 


-0.030 


-4.79466 


-4.64749 


-0.36339 


-0.035 


-4.87817 


-4.73801 


-0.36452 


-0.040 


-4.95823 


-4.82374 


-0.36696 



We use the interaction coefficients in Table H and diagonalize both Hamiltonians ( ITT]) and 
f[T2"j) utilizing the Lanczos method to determine the ground state energy at rest for the six 
considered dimers in the periodic volumes L 3 ranging from L = 6 to L = 17. We also use 
the same interaction coefficients and diagonalization method to find the ground state energy 
of the fermion-dimer systems. These energies are summarized in the tables in Appendix A. 

Now we turn our attention to the fermion-dimer scattering in a periodic cube. It is known 
that there are exponentially small corrections to the scattering energy of the fermion-dimer 
system at finite volume due to range effects. We can remove this error by extrapolation to 
infinite volume. However there is another error which is independent of the fermion-dimer 
scattering process, namely the finite volume error in the dimer binding energy. One might 
think that this error can be removed by substracting the dimer binding energy from the 
total energy of the fermion-dimer scattering system. But this is not quite correct since we 
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calculate the scattering process in the center-of-mass frame and therefore the dimer has 
some recoil momentum. The corrections to the dimer binding energy in the moving frame 
differ from its rest frame due to the topological phases for the moving dimer in the finite 
volume 0,|23|. 



III. BOUND STATE IN A MOVING FRAME 

Liischer derived the finite-volume corrections to the binding energy of two-body bound 



states for interactions with finite range [241 ] . The shift in the energy of a bound state in a 
periodic cube at rest is given by 

AE S (L) - E / d'r^Vir^if+nL), (21) 

|n|=l 

where 0oo is the infinite-volume wavefunction as a function of the relative separation f 
and V(f) is the interaction potential. Using a Galilean transformation we can find the 
wavefunction of the bound state in a periodic cube of length L moving with momentum 
2nk/L for integer k. This wavefunction in a periodic cube has a phase dependence which 
can be factorized out, 

Mr + nL) = e- 2i ™ lH Mr), (22) 

where a = m-j- / (m-j- + m±) for the general case of unequal masses and n is an integer. Each 
phase twist at the boundaries induces a shift in the binding energy of the dimer. Using 
Eqs. (|2~T!) and (1221) for the S"-wave bound state, the finite volume correction in a moving 
frame is 

AE^L)^r(k,a) AE 5 (L), (23) 

where 

1 3 

r(k, a) — - cos(2-Kaki), (24) 
3 i=i 

and AE^(L) and AE^(L) represent finite volume correction to the binding energy of the 
dimer in the moving and rest frame, respectively. 

These corrections have a universal dependence on momentum determined by the number 
and mass of the constituents. In asymptotically large volumes the corrections are exponen- 
tially small and can be neglected. But if the volume is relatively small, this shift can be 
comparable to that of the scattering process of the fermion and dimer. 

In order to calculate finite- volume corrections due to the binding energy of the dimer in 
the total scattering energy of the fermion-dimer system, let E^ d (p, L) be the total scattering 
energy with radial momentum p and Ei(L) the finite volume energy due to binding for 

the dimer with momentum 2nk/L in a periodic cube of length L. In the asymptotic 
limit L — > 00, with p scaling as 1/L, we can neglect the mixing with higher-order singular 
solutions to the Helmholtz equation. For S-wave scattering of a dimer and a fermion with 
radial momentum p and separation r between the center of mass of the dimer and the 
fermion, the position-space scattering wavefunction is 

2iivk -r 

{r\^ p ) = cY\ ^ (25) 

X 1 ' ^ {27Tk/L) 2 -p 2 
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TABLE II: Toplogical factor for six considered dimers. The upper part of table calculated by 
using the ground state energies corresponding to H\ and the lower part by using the ground state 
energies corresponding to H2. 

L 17 16 15 14 13 12 11 10 9 8 7 6 

t] 5 0.84351 0.81836 0.78862 0.75386 0.71392 0.66915 0.62052 0.56968 - 

r|-° 0.89079 0.87205 0.84912 0.82108 0.78709 0.74648 0.69916 0.64592 0.58872 

-rp 0.91790 0.90351 0.88565 0.86323 0.83517 0.80025 0.75743 0.70632 0.64778 0.58452 

t| ° 0.93511 0.92349 0.90912 0.89090 0.86769 0.83803 0.80038 0.75340 0.69669 0.63176 0.56267 



0.94645 0.93706 0.92520 0.91009 0.89060 0.86532 0.83248 0.79022 0.73731 0.67442 0.60554 



-rf 0.95466 0.94675 0.93676 0.92396 0.90735 0.88560 0.85681 0.81889 0.76973 0.70858 0.63804 0.56535 

rj- 5 0.82960 0.80247 0.77064 0.73384 0.69221 0.64649 0.59810 0.54894 - 

t| ° 0.88030 0.85976 0.83471 0.80429 0.76775 0.72470 0.67555 0.62176 0.56594 

t| 5 0.90967 0.89377 0.87401 0.84938 0.81873 0.78093 0.73524 0.68189 0.62272 0.56130 

r|'° 0.92804 0.91529 0.89930 0.87907 0.85338 0.82077 0.77979 0.72952 0.67045 0.60542 0.53965 

-rf' 5 0.94077 0.93030 0.91709 0.90025 0.87859 0.85060 0.81446 0.76847 0.71181 0.64593 0.57567 

0.94988 0.94107 0.92992 0.91564 0.89714 0.87294 0.84112 0.79953 0.74636 0.68156 0.60882 0.53618 

with some normalization constant c. The total energy Ef d (p, L) is given by 

E ^ P , L) J^M ' E I±3W (26) 
(*p|* P ) % (* 2 ->;) 2 

where M = ^2j:{k 2 — t])~ 2 . The finite- volume correction due to the binding energy of the 
dimer in the scattering process is 

Ef\p, L) - E'*(p, 00) = r d ( v ) AE$(L), (27) 

where the topological volume factor for m-j- = = m is given by 

tm = 1 y j<hiiL . (28) 

We find Td{rj) iteratively. We determine p 2 and 77 by subtracting the ground state energy of 
the dimer in the rest frame from the total energy. Using this r) and Eq. ( 128]) we find 7^(77) 
and use it to modify the binding energy of the dimer. We repeat this process until rj does 
not change anymore. For very large box length (L = 2000) we achieve a fixed point after 
three iterations. The ratio of the finite volume corrections in the moving and rest frame 
are collected in Table [III We note that this ratio is significantly smaller than one at small 
volumes. For large volumes this ratio can be neglected. The binding energies of the dimers 
in the boosted frame are summarized in Appendix A. We emphasize that the finite-volume 
correction in Eq. (127)1 has nothing to do with the interaction between dimer and fermion 
and should therefore be subtracted from the total energy before using Luscher's scattering 
relation. This subtraction reduces systematic errors in lattice calculations. We note that in 
the total scattering energy there are also corrections corresponding to the scattering process 
which we will remove by extrapolation to the infinite volume and to the continuum limit. 



8 



-1 

> -1.2 

s 

W -1.4 
-1.6 



I 1 I 1 I 1 I 1 

♦ E, = -1.5 MeV 



* * * * 



10 12 14 16 



-1.8 
-2 
-2.2 
-2.4 
-2.6 



I 1 I 1 I 1 I 



E . = -2.5 MeV 

d 



* * 



****** 



I i I i I i I i I 



8 10 12 14 16 



hi i I 1 I 1 I 1 I 1 I 



-2.7 
-3 
-3.3 
-3.6 
-3.9 



E = -3.5 MeV 

d 



"I i I i I i I i I i I 



6 8 10 12 14 16 



-1.4 



-1.6 - 



> 

£-1.8 
W 

-2 
-2.2 



I 1 I 



I 1 I 1 
-2.0 MeV 



******- 



I i I i I i I 



10 12 14 16 

L/a ,a« 



-2.1 
-2.4 
-2.7 

-3 - 
-3.3 



I 1 I 1 I 1 I 1 I 1 
E, = -3.0 MeV 



• ♦ ♦ ♦ 



I i I i I i I i I 



8 10 12 14 16 

latt 



-2.8 
-3.2 
-3.6 
-4 
-4.4 



I 1 I 1 I 1 I 



1 I 1 I 

4.0 MeV 



t i I i I i I i I i I r 



6 8 10 12 14 16 



FIG. 2: Plots of ground state energies corresponding to the lattice Hamiltonian Hi versus L. 
Circle and star represent the ground state energies of the dimers in the rest- and boosted-frame 
respectively and diamond stands for ground state energies of the fermion-dimer system. 



IV. FERMION-DIMER SCATTERING 
A. Lattice Calculation 

In order to find the radial momentum, p, in the fermion-dimer systems, we subtract the 
binding energies of the dimers in the moving frame from the total scattering energies of the 
fermion-dimer systems. The ground state energies of the six considered dimers in the rest 
and moving frame and the ground state energies of the fermion-dimer system for both lattice 
Hamiltonians, Hi and H 2 , are shown in Figs. |2]and[3l As discussed above, the difference 
between the ground state energies of the dimers in the rest and moving frame is bigger in 
small volumes. In this case it is comparable to the corrections due to the scattering process. 

The difference AE(L) is the kinetic energy of the fermion-dimer system. We use this 
to find the radial momentum p. In a naive calculation, we take 7^(77) equal to one and 
thus implicitly assume the corrections to the binding energy of the dimer the in rest and 
moving frame are equal. We subtract the binding energy of the dimer in the rest frame 
from the total energy of the fermion-dimer system in order to find AE(L) and eliminate the 
correction due to the dimer binding energy in finite volume, 

AE nmve (L) = AE fd (L) - AEi(L). (29) 

In the full calculation, we subtract the binding energy of the dimer in the boosted frame in 
order to determine AE(L) and eliminate the finite volume correction corresponding to the 
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FIG. 3: Plots of ground state energies corresponding to the lattice Hamiltonian H2 versus L. 
Circle and star represent the ground state energies of the dimers in the rest- and boosted-frame 
respectively and diamond stands for ground state energies of the fermion-dimer system. 



binding energy of the dimer, 

AE full (L) = AE fd (L) - AEi(L) 

= AEf d (L)-r d {r))AE${L). (30) 

After subtracting the corrections corresponding to the binding energy of the dimer from the 
total energy we use the Luscher's formula with this energy and calculate p cot 5q(p) for six 
different lattice spacings. We also calculate p cot S (p) for the case in which we subtracted 
only the binding energy of the dimer at rest frame from the total energy. This simply 
means we replace 7^(77) by 1. In Fig. HJ the results are plotted versus p 2 . To extrapolate to 
the infinite volume we fit a polynomial of second order to the data points. We write this 
results as dimensionless combinations multiplied by powers of the dimer binding momentum 
k. By comparing the naive calculation plots in Fig. H] (a) and (c) with the full calculation 
plots (b) and (d), we clearly see the effect of the topological phase factor. This correction 
is quite large for scattering in smaller volumes. The change in slope in plot (d) compared 
to plot (b) is expected since we tuned the effective range of interaction to zero for H2. 
From these results, we determine the low-energy parameters for fermion-dimer scattering 
and extrapolate to the continuum limit. 
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FIG. 4: Plots of pcot5o(p) versus p 2 . (a) Naive calculation for H\. (b) Full calculation for H%. (c) 
Naive calculation for H2. (d) Full calculation for H2. • , ■, ♦, ▼, x and * represent data points 
corresponding to the bound states of energies —1.5 MeV, —2.0 MeV, —2.5 MeV, —3.0 MeV, 
—3.5 MeV and —4.0 MeV respectively. 



B. STM Equation 



The STM equation for the S-wave fermion-dimer scattering amplitude Tfd{k,p;E) can 
be written as [18|, |l9 



T fd {k,p; E) 



^ M (k -E) 2 r d 1? M i<l : P> E )Tfd{k,q;E) 
3 ,Pl ttJo — k + J3q 2 /A - mE - ie 



(31) 



where k and p are the incoming and outgoing momenta of the fermion and dimer in the 
center-of-mass frame, k is the binding momentum of the dimer, and 



E = - P --- 

4 m m 



is the total energy. The inhomogeneous term 

M(k,p; E) 



1 / p 2 + pk + k 2 — mE 
In 



2pk \p 2 — pk + k 2 — mE 



(32) 



(33) 
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FIG. 5: Left: Lattice results and continuum extrapolation with error estimates for the fermion- 
dimer scattering length. Right: Lattice results and continuum extrapolation with error estimates 
for the fermion-dimer effective range. For comparison we show the continuum result obtain via the 
Skorniakov-Ter-Martirosian equation. The error bands show statistical errors corresponding to the 
fitting procedure shown in Fig. 0] 



is given by the S-wave projected one-fermion exchange. The fermion-dimer scattering phase 
shifts are obtained by evaluating Eq. ( 131]) at the on-shell point: 



Tf d (p,p; E) 



pcot 6q(p) — ip 



(34) 



By discretizing the momenta p and k the STM equation ( 1311) can be transformed into a 
matrix equation which can be solved numerically. The numerical errors in the solution of 
this equation are negligible. The effective range parameters are extracted from the scattering 
amplitude by fitting a polynomial in p 2 to 1/T} d (p,p; E)+ip. Errors are estimated by varying 
the degree of the fitted polynomial. 



V. COMPARARISON OF RESULTS AND DISCUSSION 

From the lattice calculation of the phase shifts in Fig. H]we can extract the effective range 
parameters. Our results for the scattering length, a,f d , and the effective range parameter, 
Tfd-, are shown in Fig. We analyze only the plots in Fig. H] (b) and (d) which contain 
the full calculations corresponding to H\ and H 2 , respectively. By fitting a polynomial of 
second order to each set of data we find a scattering length and a effective range in infinite 
volume for both lattice Hamiltonians. These data points are plotted in Fig In order to 
extrapolate to the continuum limit ai aU — > 0, we use a linear function. The results for the 
low-energy parameters that we get for these two independent representations of the lattice 
Hamiltonians are 

na fd = 1.162(13), nr fd = -0.041(16) for H x (35) 
Ka fd = 1.181(7), Kr fd = -0.016(16) for H 2 . (36) 

To extrapolate to the continuum limit in the lattice Hamiltonian calculations we used 
only the data points corresponding to the four smallest lattice spacings. For the other data 
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points, the Compton wavelength of the bound state is comparable to the lattice spacing. 
We estimate the systematic errors in the continuum extrapolation of the fermion-dimer 
scattering length and effective range by extrapolation to the continuum limit using only the 
first two data points and taking the interval between these extrapolation values and the 
central values obtained using all four data points as the systematic errors. The agreement 
between these two independent calculations is consistent with our estimate of the systematic 
errors. As we see from Fig. [5] the inclusion of the topological volume factor 7^(77) improves 
the accuracy, especially in the calculation of the effective range parameter. With a very 
conservative estimation of the systematic error we are able to say that the value of the 
the fermion-dimer scattering length in units of the dimer binding momentum is in between 
1.149 and 1.188. The value of fermion-dimer effective range in units of the dimer binding 
momentum is between zero and —0.057. Our final result is given by the weighted averages 
of the values in Eq. (1551) and (J3"5|) : 

na fd = 1.174(9), Kr fd = -0.029(13). (37) 

In calculating the average, we assumed that the statistical probability distribution of the 
measured variables are Gaussian and independent of each other. Using standard error 
propagation, we find the uncertainty in the average values. Our results ( 157)) are in excellent 
agreement with the continuum calculation using the STM integral equation (15T)) : 

na fd = 1.17907(1), nr fd = -0.0383(3). (38) 



VI. OUTLOOK 

We have presented benchmark calculations for fermion-dimer scattering using a Hamilto- 
nian lattice formalism and in the continuum using the STM integral equation. We obtain 
excellent agreement between both approaches. The finite- volume lattice methods presented 
here can be applied to ab initio calculations for elastic scattering of nuclei, cold atoms, and 
hadronic molecules. Of particular interest in nuclear physics are calculations of the low- 
energy scattering of neutrons upon nuclei. Soft neutron scattering upon nuclei is relevant 
to the design of container materials used in ultracold neutron experiments. They are also 
an important probe of the properties of nuclei surrounded by a dilute superfluid neutron 
gas. It is widely believed that this physical situation with nuclei in a neutron gas is realized 
in the inner crust of neutron stars. 

Calculations of deuteron-neutron scattering in the spin-quartet channel should yield sim- 
ilar results to the idealized zero-range limit presented here. There will, however, be small 
corrections due to the range of the interactions as well as spin- dependent forces. Our meth- 
ods can also be applied to neutron-deuteron scattering in the spin-doublet channel. Here 



there is an interesting connection with Efimov trimer physics in a finite volume [28M30 
As a first step in this direction, boson-dimer scattering for three identical bosons in the 
zero-range limit is currently being investigated (3lj . 

There are also some very useful applications of the topological volume factor for binding 
energy calculations. By choosing different values of the center of mass motion, one can make 
the sum over topological volume factors vanish. This can be used to remove the leading and 
even subleading finite-volume corrections to the binding energy of two-body bound states. 
This technique is being pursued in a number of recent lattice QCD studies for the deuteron 



and other dibaryonic systems [27], |32 
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Appendix A: Ground State Energies 



17 


16 


15 


14 


13 


12 


11 


10 


9 


8 


7 


6 


1.5051 


-1.5080 


-1.5126 


-1.5200 


-1.5316 


-1.5504 


-1.5804 


-1.6287 










2.0022 


-2.0036 


-2.0060 


-2.0101 


-2.0170 


-2.0287 


-2.0489 


-2.0838 


-2.1439 








2.5010 


-2.5018 


-2.5031 


-2.50550 


-2.5097 


-2.5173 


-2.5311 


-2.5562 


-2.6024 


-2.6874 






3.0005 


-3.0009 


-3.0017 


-3.0032 


-3.0059 


-3.0109 


-3.0205 


-3.0388 


-3.0743 


-3.1434 


-3.2782 




3.5003 


-3.5005 


-3.5010 


-3.5019 


-3.5037 


-3.5071 


-3.5139 


-3.5275 


-3.5550 


-3.6112 


-3.7266 




4.0002 


-4.0003 


-4.0006 


-4.0012 


-4.0024 


-4.0048 


-4.0097 


-4.0200 


-4.0415 


-4.0874 


-4.1861 


-4.3994 



TABLE III: Ground state energy of the dimers in the rest frame calculated by using H\. 

L 

El- 5 
o 

E 2.0 


£2.5 


E 3.0 


E 3.5 




TABLE IV: Ground state energy of the dimers in the rest frame calculated by using H^. 

L 

o 

E 2.0 




E 3.0 


E 3.5 


E±° 
Q 



17 


16 


15 


14 


13 


12 


11 


10 


9 


8 


7 


6 


1.5055 


-1.5087 


-1.5136 


-1.5215 


-1.5341 


-1.5543 


-1.5865 


-1.6382 










2.0024 


-2.0040 


-2.0066 


-2.0110 


-2.0185 


-2.0313 


-2.0532 


-2.0908 


-2.1556 








2.5011 


-2.5020 


-2.5034 


-2.5060 


-2.5106 


-2.5189 


-2.5339 


-2.5613 


-2.6115 


-2.7035 






3.0006 


-3.0011 


-3.0020 


-3.0035 


-3.0065 


-3.0121 


-3.0226 


-3.0428 


-3.0817 


-3.1573 


-3.3039 




3.5003 


-3.5006 


-3.5011 


-3.5021 


-3.5041 


-3.5079 


-3.5154 


-3.5304 


-3.5608 


-3.6225 


-3.7487 




■4.0004 


-4.0005 


-4.0009 


-4.0015 


-4.0028 


-4.0055 


-4.0110 


-4.0223 


-4.0462 


-4.0968 


-4.2052 


-4.4384 



TABLE V: Ground state energy of fermion-dimer systems calculated by using Hi. 



L 


17 


16 


15 


14 


13 


12 


11 


10 


9 


8 


7 


6 


pfd 
^1.5 


-1.4029 


-1.3817 


-1.3545 


-1.3189 


-1.2718 


-1.2084 


-1.1216 


-1.0003 










a 2.0 


-1.9199 


-1.9022 


-1.8791 


-1.8487 


-1.8081 


-1.7529 


-1.6765 


-1.5684 


-1.4118 








pfd 
^2.5 


-2.4313 


-2.4159 


-2.3959 


-2.3693 


-2.3335 


-2.2845 


-2.2161 


-2.1184 


-1.9754 


-1.7601 






w fd 


-2.9395 


-2.9258 


-2.9080 


-2.8843 


-2.8523 


-2.8081 


-2.7461 


-2.6568 


-2.5250 


-2.3246 


-2.0095 




pfd 


-3.4455 


-3.4332 


-3.4172 


-3.3958 


-3.3667 


-3.3265 


-3.2696 


-3.1873 


-3.0649 


-2.8773 


-2.5796 




E fd 


-3.9502 


-3.9390 


-3.9244 


-3.9048 


-3.8782 


-3.8412 


-3.7886 


-3.7122 


-3.5978 


-3.4212 


-3.1388 


-2.6686 
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TABLE VI: Ground state energies of fermion-dimer systems calculated by using Hi. 



L 17 16 15 14 13 12 11 10 9 

E(? s -1.3967 -1.3741 -1.3449 -1.3067 -1.2561 -1.1878 -1.0941 -0.9629 
E* d -1.9148 -1.8958 -1.8710 -1.8384 -1.7947 -1.7353 -1.6528 -1.5359 -1.3663 



4 d 5 

?fd 
J 3.0 
7 }d 

pfd 
^4.0 



-2.4267 -2.4103 -2.3887 -2.3602 -2.3217 -2.2689 -2.1951 -2.0895 -1.9346 -1.7011 

-2.9353 -2.9207 -2.9015 -2.8760 -2.8414 -2.7938 -2.7266 -2.6299 -2.4868 -2.2691 -1.9264 

-3.4418 -3.4286 -3.4114 -3.3883 -3.3570 -3.3136 -3.2522 -3.1631 -3.0304 -2.8269 -2.5039 

-3.9470 -3.9351 -3.9193 -3.8983 -3.8696 -3.8297 -3.7730 -3.6903 -3.5665 -3.3753 -3.0694 -2.5611 



TABLE VII: Ground state energy of dimers in the boosted frame calculated by using H\. 

L 17 16 15 14 13 12 11 10 9 8 7 6 

El- 5 -1.5043 -1.5066 -1.5100 -1.5150 -1.5226 -1.5338 -1.5502 -1.5741 - 

k 

E 2 .- -2.0020 -2.0032 -2.0051 -2.0083 -2.0134 -2.0215 -2.0343 -2.0543 -2.0855 

k 

El 5 -2.5009 -2.5016 -2.5028 -2.5047 -2.5081 -2.5138 -2.5235 -2.5398 -2.5666 -2.6106 

k 

El-° -3.0005 -3.0009 -3.0016 -3.0028 -3.0051 -3.0091 -3.0164 -3.0293 -3.0519 -3.0911 -3.1586 

k 

El- 5 -3.5003 -3.5005 -3.5009 -3.5017 -3.5033 -3.5062 -3.5116 -3.5218 -3.5406 -3.5750 -3.6372 

k 

Ei-° -4.0001 -4.0003 -4.0006 -4.0011 -4.0022 -4.0042 -4.0083 -4.0164 -4.0319 -4.0619 -4.1187 -4.2258 
k . 

TABLE VIII: Ground state energy of dimers in the boosted frame calculated by using Hi. 

L 17 16 15 14 13 12 11 10 9 8 7 6 

El- 5 -1.5046 -1.5070 -1.5105 -1.5158 -1.5236 -1.5351 -1.5517 -1.5758 - 

k 

El-° -2.0022 -2.0034 -2.0055 -2.0089 -2.0142 -2.0227 -2.0359 -2.0565 -2.0881 

k 

El- 5 -2.5010 -2.5018 -2.5030 -2.5051 -2.5087 -2.5148 -2.5250 -2.5418 -2.5694 -2.6142 

k 

El-° -3.0006 -3.0010 -3.0018 -3.0031 -3.0055 -3.0099 -3.0176 -3.0312 -3.0548 -3.0952 -3.1640 

k 

El- 5 -3.5003 -3.5005 -3.5010 -3.5019 -3.5036 -3.5067 -3.5125 -3.5234 -3.5608 -3.6225 -3.7487 

k 

Ei-° -4.0003 -4.0005 -4.0008 -4.0014 -4.0025 -4.0048 -4.0092 -4.0178 -4.0344 -4.0660 -4.1249 -4.2350 
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